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ABSTRACT 


A finite element code is developed for analysis and 
design of three dimensional truss and frame structures. 
structures are designed for minimum weight subject to 
constraints on: member stresses, Euler buckling, shell 
buckling, joint displacements and system natural fre- 
quencies. Structures are optimized with respect to member 
Seeze amd structure configuration. 

The finite element code may be used for analysis only, 
or may be coupled to an optimizer of the user's choice. 
The displacement method is used for static analysis, and 
structure natural frequencies are calculated via the sub- 
Space iteration method. 

PlowToLllegmacmpreyided for expansion to other than 
tubular frame elements, and provisions are made for the 
Peeure s2rOwr tO panel and other types of structural 
Strememus . 

Documentation is provided to facilitate use of *he 
code. A User's manual is presented with examples and re- 
Pero ieee banaclon of how this code may be coupled to 


an optimizer is also provided. 
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i. <NMRO DUCTION 


tne task of the Engineer is to provide the best solttion 
eememe  CUSGcome>'’s problem. UsWally the "best" solution is 
Ene ene that dees the job with an adequate margin of safety, 
Peeecsenetically pleasing, and is economically feasible. 
The solution may be reached through various means, but effi- 
Sememenuce of design tools, as well as efficient use of 


S important since beth add to the overall cost of 


fH. 


materials 
the product. Design optimization is one method that can be 
Meca tO Memimaize the efficieney of a structure by minimizing 
meomwement and, presumably, cost. 

Optimization of structures has had continuing changes 
Since its development in the early 1960's with an active 
area of research being elastic structures. The main goal is 
e@mdeetan StPuctural systems that efficiently perform speci- 
sled pweposes. Since most physical structures can Be 
meeelec@ Dy the »Finite element method, a computer progran 
Can be written to perform the necessary calculations to 
solve the probiem. 

The purpese of this research was to develop a finite 
element code that could be used to analyze a combined truss/ 
ereme Stmmeture and could be easily coupled to an optimizer; 


mms preveding a useful tool for desigring such structures 











as, for example, ships' masts. This code expands the pre- 
vious work by Fitzgerald [Ref. 1] on truss structures to 
the more general six degree-of~freedom truss/frame case. 

The cesign problem considered in this study is the opti- 
mization of three~dimensional statically indeterminant com- 
bined truss-frame structures under varying load conditions 
Aieeesubyect GO a variety of constraints. The objective is 
to minimize the weight of the structure where the design 
variables are member sizes and joint coordinates. Constraints 
include: maximum normal stress; maximum joint displacements; 
minimum structure natural frequencies; Euler buckling; and 
in the case of the tubular frame elements, local or shell 
muekling. 

In the present code, all gradient information is calcu- 
lated by the finite difference method. Modification of the 
code to permit calculation of gradients analytically has 
been identified as a necessary future extension. 

This document describes the use and capabilities of the 
MmEeewelement computer code to be coupled to an optimizer. 
The user's manual presented in Chapter V contains a simple 
design example in which the program is coupled to the CONMIN 
optimization code [Ref. 2]. Additionally, guidelines for 
coupling the code to an optimizer of the user's choice are 


presented. 


RZ 





several examples demonstrating the program under a 
variety of conditions are presented. Conclusions and reconm- 


mendations for future work are given. 


| 
CW 
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A. INTRODUCTION 

When the finite element method of analysis is used to 
design optimization, two objectives must be kept in mind. 
First, the number of analyses for the structure should be 
kept to a minimum. Second, the amount of gradient informa- 
tion required during the design process should be minimized 
to shorten run times and minimize computer storage 


requirements. 


B. STATIC ANALYSIS 

Initial formulation of the problem must include approxi- 
mate member areas in the case of truss elements, and for 
frame elements, characteristic dimensions (for tubular members: 
mean diameter and wall thickness); material properties (which 
may be different for each member); a set or sets of external 
loads; any non-structural attached masses; and specified 
jetneesupport conditions. 

The analysis for the stresses and deflections at the 
Jedieswmust satisfy the conditions of equilibrium of forces 
at the nodes and geometric conditions of compatibility of 
deformation. In this analysis the structure is assumed to 
behave in a linearly elastic fashion. The weight of an 


individual member is not inherently included as part of the 


14 





egeecttied toad conditions; but, as an option, half of the 
weight of each member may be applied at the member's end- 
points as an additional load in the negative Y-direction. 

For this analysis, the following assumptions are made: 
truss and frame members are treated as discrete entities; 
truss elements have three translational degrees of freedom 
at each node and are treated as pin-connected; frame elements 
have three translational and three rotational degrees of 
freedom at each node and are treated as fixed-fixed beams; 
and loads and reactions are applied at the joints as shown 
Miesicgure 2.1. 

The Displacement (Stiffness) method for finite element 
analysis [Ref. 3] is utilized where 


C34) 


RA 


Ue 


4 


and where is the global stiffness matrix, F is the vector 


wx 


Bmevectors Of applied loads, and u is the vector or vectors 
of displacements. The method used herein is an extension of 
that described by Felix and Vanderplaats [Ref. 4]. By 
applying the constitutive stress-displacement relationships, 


stresses in the elements may be recovered. 


C. DYNAMIC ANALYSIS 
When constraints on the system's natural frequencies 
Mees to be considered, the design process requires the solu- 


tion of an eigenvalue problem. The sub-space iteration 


aS 
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meoieem2.! FORCE ORIENTATION CONVENTION FOR AN ELEMENT 
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method of Bathe and Wilson [Ref. 5] is used to solve for the 
desired number of lowest eigenvalues and the associated 
elgenvectors. This method is reasonably efficient for a 
small number of lowest frequencies for a large problem, and 
is well suited for re-analysis when small changes are made 


in the design. 


Pee GRADIENTS 

Gradients are currently calculated with respect to member 
sizes and coordinates only by finite difference techniques. 
Inclusion of the capability to calculate gradients analytic- 


ally is identified as a needed extension to this work. 


ie? 





it. OPE cA teen 


A. INTRODUCTION 

Mies primery Objective of structural optimization, is to 
design systems that efficiently perform specified purposes. 
Selection of a specific optimizing algorithm must take into 
account the following: 1) the structure should be analyzed 
as few times as possible; 2) the algorithm should minimize 
the amount of gradient information required; 3) the algorithm 
should provide reasonable assurance that an optimum or near- 


optimum design will be reached. 


Be GeNeRAL FORMULATION 
The general statement of inequality constrained minimiza- 
tion is as follows: 
Minimize 
EC® (4.1) 
Sie gect to: 


S.(¥) < 0 j=ljm (4.2) 


ei 2 ae be T=i5n (4.3) 
i-°1il- “1 


where F(X) is the objective function to be minimized. The 
mUNIGT TORS SNS are the set of inequality constraints to 


be met. The vector X contains the design variables. 


dls 





The inequality constraints, 63¢%) < 0 F221 5m) oniet be 

eeetoried for the design to be accepted as feasible. Side 
; = & u 

Seomsmmaznts , Xe and xX; 5 are lower and upper bounds on the 
design variables. The objective function must be minimized 
as much as possible while still satisfying all inequality 
S@eeerain=s. in the case that it is not possible to satisfy 
all constraints, the most nearly feasible solution must be 


Found. Felix and Vanderplaats [Ref. 4] is an excellent 


Bemres fOr the basic structural design formulation. 


ee DESIGN VARIABLES 
imemveetor X contains the design variables; in this case, 
member characteristic dimensions for frame elements, cross- 
Sectional areas for truss elements, and spacial coordinates 
Pewee joints. the user may elect to optimize the structure 
weight with respect to any combination of the above variables. 
For truss elements, where normal stress 1s dependent only 


an element 


fu 
O 
tty 


Meogmuenecemesmitude of the cross-sectional are 


ct 


TLeLence. © 


t-> 


ty 
4; 


mae notetie cistribution of tha Seo 


area, it 
use the area as the design variable as long as the Euler 
buckling stress can be related to the cross-sectional area. 
For frame elements where stresses are dependent on area 
@istrzomtion, meme flexibility is allowed by varying two 
P@iserersstic dimensions and allowing the code to calculate 


mame Sectiem properties from these. 





In the case of tubular elements, the characteristic 
dimensions are mean diameter and wall thickness; from which 
area, bending and polar moments, and maximum radial dimen- 
Sions are generated. 

The optimum geometry problem requires that the joint 
coordinates be design variables. The x, y, and z coordinates 
of a joint are treated as separate design variables. 

In many cases it is desirable to link a set of design 
variable together to preserve symmetry, limit the number of 
variables to be solved, or to limit the number of unique 
elements to be manufactured. The code has provisions for 
design variable linking by which two or more variables may 


be linked in equality or some fixed ratio. 


Pee OBeECTIVE FUNCTION 
Mie op jective function under consideration is weight 
N 
F(X) = Ee ioe, (4.4) 
Rew 
where 9; is the weight density (in consistent units) of the 
material of the ith element, A. is the cross-section area 
Seeene ath element, and Ls is the length of the ith element 


and NE is the number of elements in the structure. 


Ecos’ RAINTS 
This code is designed to accommodate constraints on maxl- 


mum normal stress, Euler buckling, local or shell Due wl ime. 
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maximum joint displacements, and minimum natural frequencies 


Or the structure. All constraint values are normalized. 
i = Semess: 


Yer eeiss elements, maximum normal stress a, is 


Frame elements; tensile stress O, ,» maximum dendin 
Shear stresses On and oO, are calculated. Maxithum and 
Tenomum norma! stresses are 


G = max magnitude of o, + ¥(o vet On es 
max al 





fins 9) 
Q 
+ 
oN 
Q 
+ 
f Bey 
Q 
J 
Nr 
cn 
[= 
O> 
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m= Min magnaizude o 
min 


NO 


ct 
ty 
O 
DP) 


The upper and lower bounds on stress may be dirreren 


each member, Dut are taken to be the same for every loading 


comcdition. 


ams pre = 5 pean ae. . sans * > 7 7 _ 
Phe Samess at which local or thin sneit, SsucxLang 
Seames 15 given ey: 
OC - a 0.456. oe jy) 
lus or as 
2 i 
M1 | 7 a4 
Mere the subscript i corresponds to the member numter, £E. 
is Young's modulus, D_ is the element's mean Clameter, and 


7 
—7 


me tie efement's wall thickness. 


ct 


[- 


oe) 





F. DESIGN VARIABLE BOUNDS 


ry 


Side constraints are imposed on the design variables 
es: 
& u 


ores D. < CD; (4.8) 


: ; : ; . Y u 
Gees che characteristic dimension and CD. and CD. are 
mee mMinamum and maximum aellowable characteristic dimensions 


of the ith element, and are taken to be the same for all 


load cases. If, as is the case in a tubular element, geo- 
metry dictates that some relationship between the design 
Variables cannot be exceeded, i.e., thickness cannot exceec 


mrommeam Giameter, the user must arrange the bounds to pre- 


mlie= Such an Occurrence. 


eee, COLL IMUM GEOMETRY 


Joint coordinates are treated as design variables with 


oT 


a separate design variable for each coordinate direction 


= 
( 


of a node. Coordinate design variables may also be linked. 
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ie EROGwealt PEATURES 


A. INTRODUCTION 

Computer codes each have their own features and formats 
With which the user must become familiar if he is to use 
the code easily and efficiently. The SADX code developed 
in the course of this research has been designed with 
Pomtrous Options which are explained in this chapter. 
Chapter V contains a User's Manual with sample data for a 
typical problem thet might be solved with this code: a 
truss-braced cantilever beam. This problem along with 
other numerical examples will be presented in detail with 
results in Chapter VI. 

The SADX code was written to be used as a stand alone 
analysis program or as an analysis code that could be easily 
be coupled to an optimizer (of the user's choice) through 


simple modifications to the main driver program. 


Be FEATURES 

When the user supplies member areas, section types, 
Characteristic dimensions, material data, connectivities 
and joint coordinates, along with a set of program control 
parameters; the analysis mode will calculate the weight of 
the structure. The addition of one or more sets of loading 


parameters wiil result in the calculation ©: Tesultant Jorn: 
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displacements, member stresses, and/or forces for each load 
case along with the desired number of structure natural 
frequencies and modes. Design variables may be chosen as 
truss element areas, frame element characteristic dimen- 
sions, and joint coordinates. When coupled to an optimizer, 
the code will minimize the weight of the structure and print 
the final optimization information. The code is designed to 
be as simple to use as possible while retaining the flexi- 
bility for use on a variety of problems. 

The code'’s modular construction allows the user to use 
frame element cross-sections other than tubular elements. 
This is done by reading in two characteristic dimensions for 
each frame element along with a section type identifier. 
Subroutine SADX85 is called to calculate the section proper- 
ties, area, maximum radial dimensions, and bending and polar 
moments of inertia. The user may augment this subroutine to 
calculate section properties for whatever section type he 
may choose to work with. 

Many computer codes require the definition of an auxil- 
iary node to orient the principle axis of a non-axially 
symmetric element. In this code it is assumed that the 
element's local Z-axis is parallel to the global x-z plane. 

Pseudo-dynamic storage is used to allow storage of most 


data in one integer and one real array for more efficient 


use of storage. 
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Various print control options are available to tailor 
the printed output to match the user's desires. 

Design variable linking is available to allow elements 
and joints to be grouped to maintain symmetry, limit the 
number of independent design variables, or reduce the 
variety of member sizes generated. 

Optimization may be performed with respect to member 
Size, with respect to structure geometry, or with respect to 
BOE. 

Structures may be comprised of truss elements, frame 
elements, or a mix of the two types. 

structures may be optimized for multiple load cases with 
constraints imposed upon any combination of maximum normal 
stress, maximum joint displacements and rotations, Euler 
and local Buckling, and a specified number of minimum natural 
Frequencies in free vibration. Separate displacement con- 
straints may be imposed for each load case. Both consis- 
tent and lumped mass options are coded. Either forces or 
stresses or both can be output. The user may decide whether 
to include the structure's weight and the fixed masses as 


loads applied to the structure. 


Ge EXAMPLE 
ijemrtollowing example of the truss-braced cantilever 
beam presented in Tables (I-V) demonstrates some of the 


options available in the code. 
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A. INTRODUCTION 

In developing any computer code for engineering analysis, 
it 1s necessary to additionally develop concise, easily 
understood documentation. This SADX USER'S GUIDE is written 
to be easily understood by the user having only minimal know- 
ledge of the FORTRAN language. The format follows that of 
the optimization code, COPES/CONMIN [Ref. 2]. 

This chapter is devoted to acquainting the user with the 


code and necessary input data. 


Pree estGN EXAMPLE 

The simple example of a four-element combined truss/frame 
structure is used to demonstrate some of the features of the 
BeawA program. 

Giemecructure is shown in Figure 5.1, and consists OF 
two tubular frame elements along the x-axis with two truss 
braces to the y and z axes from the joint between the frame 
elements. A non-structural fixed mass is attached at the 
cutboard end of the second frame element where two loacs are 


applied. 


C. SADX DATA 
The SADX program reads data from unit 5 and writes out- 


SUe,on unit 6. Units 30 and 40 are used as Seratcn @l1ies. 





ORIGINAL 





Figure 5.1 TRUSS-~BRACED CANTILEVER BEAM 


oe 





The scratch file numbers may be changed by changing two cards 
at the beginning of subroutine SADX0l1. The SADX program has 
the capability to read unformatted data. The following sec- 
tions define the data which is required by SADX. The data 

is segmented into "BLOCKS" for convenience. 

SADX data begins with a TITLE card and ends with a END 
card. Comment cards may be inserted anywhere in the SADX 
data stack prior to the END card, and are identified by a 
doilar sign ($) in column 1. Data may be of either the 
"T10" or "F10.0" type or may be free-format separated by 
commas with no imbedded blanks. Lines of formatted and 
unformatted data may be intermixed. 

iesormatted Data 

femicdes are of “T10" and "“Fi0.0" type. "I" formats 
Meee oe might justified, and "F" formats must have the deci- 
mal point. The number of cards read per data block is listed 
feeene bottom of each block. 

2. Unformatted Data 

While the USER'S MANUAL data sheet defines SADX data 
in formatted fields of ten, the data may actually be read in 
a Simplified fashion by seperating data by commas or one or 
more blanks. If more than one number is contained on an un- 
formatted data card, a comma must appear somewhere on the card. 
Tewespomential numbers such as 2.5+ti0 are read on an unfcr- 
matted card, there must be no imbedded blanks. Unformatted 
cards may be intermingled with formatted cards. Real numbers 
on an unformatted card must have a decimal point. 
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EXAMPLES 
Unformatted data 
Poel, Sto, =5 ,2 
Equivalent formatted data 


eo |. 10 20 30 40 90 60 70 80 
é iy S72 pies to -5 0 2 





Unformatted data 
2 

28 

263 


Equivalent formatted data 


eel 18 20 30 40 20 60 70 80 
2 
2 3 

Le) 


NOTE: The third line of data contains no commas and is 
therefore assumed to be already formatted. 

Placement of more than eight unformatted data on a 
card will create two (or more) formatted cards as required. 
Fields of zeros will be created if more data are required 


mifaneeeeeo riled on an unformatted card. 


D. CONSTRAINTS 

Constraints are calculated and stored in the G vector 
Boies in the following chart. The total number of 
eonstraints 

NGON=NFREQ+2*“NDSPLC+NLC*(2*NEB + uU“NEF + NEB + 2*NEF) 


rTreq displ stress SEress Duel buck] 


3 4 





Where NCON is the total number of constraints, NFREQ is the 
number of frequency constraints, NDSPLC is the number of 
displacement constraints, NLC is the number of load cases 
imposed, NEB is the number of bar or truss elements, and NEF 
is the number of frame elements. When any of the constraints 
are missing from the G vector, all constraints are moved up 
in the vector. For example, if there is no frequency con- 
straint, then a displacement constraint would fill the first 
location of the G vector. 

Seio@aatiNTS ARE STORED IN THE G VECTOR IN 


THE FOLLOWING ORDER: 


Maer e@eEREQUENCY CONSTRAINTS 
2*NDSPLC JOINT DISPLACEMENT CONSTRAINTS 
Meese ONSTRAINTS ARE STORED ELEMENT BY ELEMENT 
FOR A GIVEN ELEMENT CONSTRAINTS ARE STORED BY LOAD CASE 
FOR A GIVEN TRUSS ELEMENT AND LOAD CASE, 
CONSTRAINTS ARE STORED: 

(Eee STRESS LOWER LIMIT 

Metetwhe STRESS UPPER LIMIT 

EULER BUCKLING STRESS LIMIT (IF APPLICABLE) 
FOR A GIVEN FRAME ELEMENT AND LOAD CASE, 
CONSIZAINTS ARE STORED: 

NEAL STRESS AT LOW NODE LOWER LIMIT 

N@BMAL STRESS AT LOW NODE UPPER LIMIT 


Mette STRESS AT HIGH NODE LOWER LIMIT 
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NORMAL STRESS AT HIGH NODE UPPER LIMIT 
EULER BUCKLING STRESS LIMIT (IF APPLICABLE) 


LOCAL BUCKLING STRESS LIMIT (IF APPLICABLE) 


ae oxXAMPLE 

The initial configuration of the braced cantilever beam 
is shown in Figure 5.1 Stress constraints are imposed as 
well as constraints on Euler and local buckling, displace- 
ment, and first fundamental frequency. A non-structural 
mercd Mass 1S applied at the tip of the beam, and two load 
conditions (Pl, and P2) are imposed. The structure's own 
weight will be considered as an imposed load as will be the 
fixed mass. 

The linking of design variables is demonstrated by link- 
ing the mean diameters of the two frame members (Dl, and 
D2 with D1=D2). Member size variables are: truss member 
areas (Al, and A2), frame member mean diameters (Dl, and 
D2), and frame member thicknesses (Tl, and T2). Member 
Sizing DESIGN VARIABLES are: XA(C1)=Al, XA(2)2A2, XKA(3)=D1=D2, 
Peach )y=T1, XACS)=T2. 

Geometry variables are the attachment points (joints 4 
and 5) of the two truss members on the y and z axes (74,25). 
Coordinate DESIGN VARIABLES are: XC(1)=Y4 and XC(2)=2=Z5. The 
structure's weight is then optimized with respect to member 
sizes and structure geometry for a total of seven design 


Variables. 





i eeeesoperties/Conditions 


Iwo material types are used: type 1, aluminum, is 


used for the truss members; and type 2, steel, is used for 


the frame members. The weight densities of the materials 
(o) are: 

Gygee 1 = 0.1 Wey sag 

mye 2 = 0.3 ib/in? 


ive Young's moduli of the materials (E) are: 


MOM OE+S psi 


crf 
TT 


type 1 


typem2 o£ Ze 205 PC ses. 
ieew en sSeructural fixed mass attached at the tip of 
semucture weights 250 15. The applied loads are: 

me mOO0.0 Ib in the +ty direction 

Ee MOO bo In the =z direction 


The acceptable maximum normal stresses are: 


eee  -25000 psi < 


C < 25000 psi 
= ia — 


oO 


myee 2 —-36000 psi are 


[A 
LA 


36000 psi 


Displacement limits, imposed upon joint number 3 (the 


Mmemeedem load case in the direction of loading, are: 


load case 1 y-direction tia, 30> 2s 
load case 2 Z-direction if 13 See 


Poumes., placed on the positions of joints 4 and 5 al 
PamemoeZ axes, are: 
fied 4 y-coordinate 0:0 @henes: to. 200.0 


Z-coordinate O20 a3enes Fou 0e.. 6 


cm 


quejabelns 


the 


t iep) 


imehes 


Linenes 


Memon natural frequency of the structure #s constrained 


cr 


memeoe greater than 1Hz. 


cm 





2. Input Control Parameters 
The following input control parameters are given 


for ease of following the example: 


NEB=2 NEF=2 NJ=9 hCw=3 
NMT=2 LDV. CEe=2 NDJ=2 NEUBC=1 
LBUCK=1 NFREQ=1 NFMASS=1 NEIG=1 
Peel = 2 NPRI=0 NLC=2 NDSPLC=2 
NSERES =2 NSTW=1 NFMW=1 


Table VI is a listing of commonly used nomenclature. 

The following USER'S MANUAL is divided into blocks 
A through P. Appearing directly below each data field line 
are the parameters for the TRUSS-BRACED CANTILEVER BEAM 
Peemple. it is important to note that the user may choose 


any consistent system of units. 
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A 
BL 


BU 
Sank DIM, 2 
DER 


GRAV 
Ov CabC 


20.4 


teed 


ee ee 


COMMON VARIABLE NOMENCLATURE 


“Usmpon's CROSS-SECTIONAL AREA 
LOWER BOUND ON DISPLACEMENTS 
UPPER SOUND ON DISPLACEMENTS 
CHARACTERISTIC DIMENSIONS FOR FRAME MEMBERS 
POR LSECT.EQ.1: MEAN DIAMSTER €2ND THICKNESS 
es a 1=X, 2=¥, 3=Z, 4=ToOt adout x, 
adoutey , é6=rot about z. 
YOUNGS. MODULUS 
CONVERGENCE [TOLERANCE OF EIGENVALUE 
SOLUTIONS eee eee 
LOMER BGUND ON FIRST,SEGOND, ETC. 
Ss POM NATURAL FREQUENCY 
LOSD EeaCss APPLIED IN THE X DIRECIION 
Gime FORCES APPLIED IN THE Y DIRECTION 
LGMD FORCES APPLIED IN THE Z DIRECTION 
ACCELERATION DUE TO GRAVITY ee 386.4) 
THe OPTIMIZATION OPERATION IDENTIFIER 
1 FOR OPTIEN UM OSMBra SIZE ONLY 
2 FOR BOTH OPTIMUM MEMBER SIZE AND GEOMETRY 
So PURMO- TIM US GHOMETRY ONLY 
COMA C a ene co Swe Peeks Lr NON-ZERO @iHe -X=DOr 
PoeeouotrRal NED 
Gon aa mit TDENTIPIERS ITS NON-ZERO THE Y-DOF 
IS CONSTRAINED 
GONSDHAGET Dh NT PEPER. IF NON-ZERO THE Z-DOF 
DS COWS TRAT NED 
CONSE oeNY GeeNT Liege. ITF NON-ZERO THE 
Poeas LON TAS OUT THE X-AXES IS CONSTRAINED 
Colon ont TDENTIFIER. IF NON=ZERO THE 
Roma t SON ABOUT THe Y-AXDS 1S GONSTRAINED 
COS RANE PLU aN. cE LOR. Le NON-ZERO TERE 
RelnT GON 250UT THE Z-AXES IS CONSTARATW=D 
JOINT NUMBER (GLOBAL) 
Cue tem GUL En SUCKLING COSPPICTEnNT SOS. Sak 
BLEMENTS 
PoCeL SUCKLING CONSTRAINT IDEATIFiIZR 
ie UCL. ne Oe LOCAL BUCKLING CONSGRRENTS 
teooees leo bbe TO TISULAR FRAME MEMBERS 
LOAD CONDITION 
LUM2SD WaSS OPTIONS 
iF LMASS.2Q.9 CONSISTENT MASS MATRIX USED 
TeeecnaSS.Ns.0 LUBPED BRASS MRT RIX USED 
ELEVEN: SUM SER 
Cxwoce oe Oer ON TYet TUSNtiri zr? 
LSECT.20.1 INDICATES A TOBULAR MEMBER 
THE NUMBER OF CONSTRAZNED JOINTS 
THE NUMBER OF JOINTS WITH (OR LINKED TO) 
BESUMGN VARI ASLES 
TROSS MEMBER AREA DESIGN VARIABLE NUMBER 
COORDINATE DESIGN VAREABL= NUMBE2 
PRAME MEMBER CHARACTERISTIC DIMENSION 1 
DESIGN VARIABLE NUMBCR 
FRAME MEMBER CHARACTERISTIC DIMENSION 2 
DESIGN VARIABLE NUMBER 
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DATA BLOCK 8 
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NOTE: DEFINITIONS OF PROGRAM INPUT CONTROL 
Pamamere RS ARE LISTED ON NEXT PACE 
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CONTENT 


PURSE “GARD 
NEB~number of bar elements 
NEF-number of frame elements 
NJ-nember of joints 
MCcamunoeweot Conscrained joints 
NMT-number of seperate material types 
DOVELC—-d=esign variable concrol parageter 
If (IDVCLC.EQ.1) NDV=NDVAR1 
Optimizes wrt member siza only 
If (IDVCLC.EQ.2) NDV=NDVAR1+ NDVAR2 
Optimizes wrt memper sizes and g¢ometry 
If (IDVCLC.EQ.3) NDV=NDVAR2 
Optimizes wrt geometry only 
NDJ-total linked and design variabie joints 


i.e. number cf 'movable’ joints) 


SECOND CARD 
NeUoe=-Eyncr buckling Constraint 2L¢entiiis: 
TET (NeW ee. 1S. 0) LULER buckling con 


Woe) be aiposed on bar elsmer cs 


+e | 
ct 
ij 
tu 
4) 
) 
ct 
'' 
th 
i? 
(D 
| 


LBUCK-local buckling constraz 
moe VeKoNr sO peLoca 2 sb 
#ill be imposed on tubular fr 
NPRB2QO-number of frequency constraints 
NFMASS-number of fixed nasse 
Secuc-u7s 
NEIG-number of precis? eigenvalues #9 


be avaluated 
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NEIGI-number of eigenvalues to be avaluated 
DEFAULT=min. of (2*NEIG , N®IG+3) 
MPRI-print control identifier 
NPRI.,ne.0 input info not printed 
Men1.26.5 RA/TA/LOCR/LICIL atrevs 


Will be printed for debuggirg 


THIRD CARD 
NLC-number of load conditions 


NDSPLC-numbe> of displacement constraints 


NSTRES-ferce/stress print identifier 
If (NSTRES.EQ.0) stresses are printed 
MemGtcot ho S.50. |) horees ere ortrtsd 


Hi 
rh 


(NSTRES.EQ.2) both 2r2 printed 
NSTW-structure weight identifier 
If (NSTW.NEZE.0) the structure's waight is 
considered as loads 
NOMW-fixed mass weight identifier 


TE (NPMW HE.) the fixed Gaasses a 
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ESCRIPTION: Dynamic Analvsis Information 





on lil 


FIELD CIN TENTS 
1 LMASS-Lumped mass option (if LMASS.N&.0O) 
the lumped mass matrix is used. 

Z GRAV~accleration due to gravity 
(default=386.4 inches/sec?) 
3 EPSEIG-conversence tclerance on ¢igenvaius 


S Odes © Rewetdeiea ul +=.0109 1} 
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BA=- jOLn= COOh1i nate 2 


X-x cocrdinacts 
m—Y coordinat? 


Z-zZ ccordinai= 


HOTE: Number cf 
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Ome: thts block TE NDG=0 in BLOCK B 


DESCRIPTION: Coordinat? Design Variable Linking Data 





Beater 
Wa 106 1020 


IX-design variable associated with x coord. 


LT¥-design variable associated with y coord. 


, 
zZ 

3 IZ-~design variable associated with z cocri 
4 PCX-participation coefficient of x-coord. 
5 


PCY-particip2ticn coefficient of V=CoOgG. 


6 PEZ—DaL moeinde On copa. wes =n~ IO. Z~coord. 


NOTE: Number of cards r2ad=NbDJ 


Fotma~ and "Example 


=a SIGMIN ISIGMAX IKEULERIPOTSSN] FORMAT 








FIELD CONTENTS 
1 p-Young's Modulus 
2 RHO-material density 
3 SIGMIN-mininum allowable stress 


ep 
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eS 
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SIGMAX-naxinum allowaple 


j+ 


Ss 
KEULER-Zuler buckling cosfficient 


NO MN £& 


POITSSN-Poisson's ratio 


NOTE: Number of cards read=NMT 


4 8 





C8Bt 261s ‘pieck 82 Mes=0 2n BLOCK 8B 
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1 LNO~element nuaber 
2 NODE I-qlobai mumber associated with the 
S 1st eie 


4 MATCOD=-material type of this slement 
YDSG1-desian variabl2> number associated wWltn 
“his element's areca 
A-member cross-sectional area 
LSpel-csocs-Stction tyoe identifier 
ESsC?TorO.t indicates ae tubutel gemper 








GCtese ea nas block if N2F=0 ir BLOCK 3 
Meecki STION: Frame Elemen= Infozrnation Faotrmat and 
Example 





FIELD CONTENTS 
1 LNO~element nunber 
Z NODE1l-qlobal number associated with the 
element's Ist node 
3 NODE2-gqlobal aumber associated with the 


element's 2nd node 
MATCOD-material «ype of *his element 


NDSG3-design variable nunaber associated with 


the elemen*'s Ist characteristic diimensior 
6 NDSG4-design variable number associated with 
“he element's 2nd characteristic dlimension 


7 Motel -cCECes-Secticn type identifier 
moanerer Ost to d.ca tes a tubula= nemoce 


NOTE: Number of cards read=NcP 
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DATA BLOCK J 
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Omit this DECCK Te eEE=0 =o BLOCK 2. 


Format and Example 


| FORMAT? 
I10,6F10 





FIELD CONTENT 
1 WLG=numbDeT Of bead=dejo12ts £92 =his Lead 
Cond 22200 
1 JN-joint number 
2 Px 
3 mi TOrces in fn] Lah, 6 1etbeer iors 
4 ra 
5 al 
ih 


~~ ON 
ae | 


NOTE: Number 
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Number of sets of c 





DATA BLOCK kK 
Omi< this block if NFMASS=0 in BLOCK 8B 


MeocRIPTION: Fixed Mass Information 


Format and Example 
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FIELD CONTENTS 
1 JN-joint number 
2 PeeSSUOTE- dass at jein: (ia)) sa tosceuun ts 


NOTRE: Number ct cards read=NFMASS 
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DATA BLOCK L 


Oma thes block 2f£ TDVELZ=3 


DESCRIPTION: Design Variable Information 
(MEMBER SIZE VARIABLES) 


Format and Exampile 


XA (1) XA (2) FORMAT 


SF10.0 





20.0 20.0 20.0 2.0 2.0 








FORMAT 
SF 1Om0 





XAL (1) | XAL(2} 





XAU (1) XAU(2) 


(is) 16, g) a, Lae Des 4.0 


KA-initial value of area design variabls 


NO 


XAL-lower bounds on area design variables 


XaU-upper bounds on arsa design varieb ls 


NOTE: read one value of XA, XAL, XAU for ¢€a 
independent member size variable defined in 
BLOCKS G and H 


Number of cards read = as required 
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Format and Example 
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23 mey-lowes bounds on coord. design var2a5les 
3 Meg-upeper bounds on coord. design vana25.2s5 
Weme=e oad one Value of XC, XCLZXCU EOr cach 
independent coord. variable defined in BLOCK D 


Number of caris tread =as required. 
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PESeeeP TION: Joint Displacement Constraint Informa*icn 
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FIELD CONTENTS 
1 JN-joint number 
2 NIR-direction 


1=x, 2=y ,3=Z iisplacemen* 
Uae 5=7 ,O=Z = Otation 
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- Bi-lowez bound on displace 
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5 BU-upper bound on displacen 


NOTE: Number of cards read= NOSPLC 
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VI. NUMERICAL EXAMPLES 


A. INTRODUCTION 

Design of three-dimensional truss and frame structures 
ape presented herein and the corresponding numerical results 
are summarized to demonstrate the use of the SADX code. 


Biemexamples begin with the truss-braced cantiliver 


beam. 


B. EXAMPLE 1: TRUSS-BRACED CANTILIVER BEAM 

The simple truss-braced cantiliver beam, as shown in 
Peeure 6.1, has been previously used for the SADX USER's 
MANUAL example. It consists of two steel tubular frame 
members with a common outer diameter and different wall 
thicknesses arranged as a cantilever beam along the X-axis. 
There is a fixed 250 pound mass at the tip of the beam. 
Two aluminum truss members are attached from the beam mid- 
BetncecO DOInts on the Y and Z axes. This structure is 
designed for optimum member size and geometry under a set 
©6f two load conditions and subject to constraints on maximum 
stress, maximum joint displacement, Euler and local buckling, 
and minimum structure natural frequencies. The weight of 
the non-structural fixed mass and the structure's own weight 
are imposed as loads. There are five member size and two 


GGerdimate design variables, and a total of 41 constraints. 
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The number of analyses required for this design is 116 using 
2.7/5 seconds of CPU time and terminating on the 12th itera- 
tion. Of the analyses conducted, 84 were required for the 
caiculation of gradients. The weight of the structure, 


including non-structural fixed masses, is reduced from 9542 


peunds to 338 pounds. Results are given in Table VII. 


iCXAMPLE 2: TWO-TIER 3-D PORTAL FRAME WITH TRUSS X-BRACES 
Une we@e-tier three-dimensional portal frame with truss 

member diagonal braces, as shown in Figure 6.2, is a symme- 
tric moderately sized case that can be analyzed easily by a 
Vemiety Of Other codes. There ere 20 truss elements, 16 

Frame elements, and four non-structural fixed masses. The 
Ieterial used is steel. This structure is designed for 

Optimum member size and geometry under a set of three load 


Sone#t2ons and subject to consfraints on maximum stress, 


ieeermum Joint displacement, Euler and local buckling, and 
minimum structure natural freauencies. The weight of the 


non-structural fixed masses and the structure's own weigh 
@re imposed as loads. There are 10 member size and five 
S@O"m@inace Césign variables and a total of 493 constraints. 
The number of analyses required for this design is 385 using 
M7 Weeecemes Ct CPU time @nd terminating on the 21st itera- 
20m. Of tne adhalyses coneucted, 305 were required for the 
StU clon of sareadiencs. fhe weight ofthe structure, 
Wieimding non-structural fixed masses, is reduced from 9302 


Pemees co 1462 pounds. Results ere given in Table VII. 
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Sree seirhe 3: DD-963 FOREMAST 

Example three is a redesign of the forward mast on the 
DD-9653 of SPRUANCE class destroyer. The DD-963 foremast has 
Deen chosen as the third test case for the following reasons: 
meeeee structure is typical of the masts found on many com- 
Daeants in the United States and other navies, 2) high top- 
side weight reduction is desirable from a stability viewpoint 
for any ship, 3) the structural members are predominantly 
tubular, 4) the problem is sufficiently complex to make 
conventional design methods cumbersome, 5) the member size 
and loading information is available. 

Bresseructure aS shown in Figures 6.3 through 6.6 is 
constructed of 172 aluminum frame members. The material 
used in the structure is 5086 H32 aluminum, an alloy with 
Meaerabte strength, good weldability, and good corrosion 
mesistance. 

This structure is designed for optimum member s1izZe under 
Peotmete 1oad condition and subject to constraints on maxi- 
mum member stress, maximum joint displacement, Euler and 
Meeal buckling. 

Some structural simplifications are made. The weights 
of mast~mounted radars, antennas, and platforms are modeled 
by a series of fixed masses which are imposed as loads along 
Mery the structure's own weight. The forces due to halyareds 
and wire antennas are applied as loads. Inertial forces due 


Bemeniips motion are calculated for the initial design point 
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Figure 6.5 DD-263 FOREMAST FORWARD SIDE 
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Vil. CONCLUSIONS AND RECOMMENDATIONS 


oe CONCLUSIONS 

An existing finite element code was expanded to encom- 
pass the more general case of frame members; i.e., six 
degrees of freedom per joint. Combined truss and frame 
eemuetures were designed for minimum weight with multiple 
load conditions considered. 

The displacement method was used for static analysis 
anc the subspace iteration method was used for eigenvalues. 
several examples were considered. In every case the 

code worked as an analysis tool, and significant weight 
reductions were obtained with the coupled optimizer CONMIN. 
The SADX design code has been shown to be a useful tool 


fOr ship mast optimum design. 


pee Sa COMMENDAT TONS 

Mies: ollOwing recommendations may be of value for fuvure 
monk . 

1. The routines necessary to calculate gradients analyt- 
ieally should be acded to the code. 

Queene, coce should be extended to include plate and 
BPeeareelements and a library of member cross-sections. 


See out of core equation solver should be acded. 


~~ 
Ca) 





* 


eerie adility to specifyemilwipliems for applyigs mer- 
Gial loads along the three coordinate axes should be added 
Meer asiiOn similar to that used for appdying stmmetume's 
Own weight és loads. Such an addition would simolify dynami 
load analysis ane Ces tan 

9. The method of gradient calculation should be deven- 


‘ 
_ 


Sememem specific sradients requimed [Ref. 7] and lRef. €]. 


6. Gradients of frequency constraints would benefit from 


Minmeeacde still exists. 
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APPENDIX A 


Dyes wees Linas 


INTRODUCTION 
This appendix contains the data files used to create 
test cases in Chapter VI. Additionally the data file 


the USER'S guide in complete form is presented. 





TABLE X 


DATA FILE TRUSS-BRACED CANTILEVER BEAM 


$RSRHK AHH] AY AS AKA? HD ERS OH AY RRS HH GH ERE ARES HEAR ERAES, ARKAEHKERET EKEEE ESN 


PoLeck A. TITte 1 CARO FORMAT H 
TRUSS“BRACED CANTILEVER BEAM: EXAMPLE 1 a 


$R SRS IGS] HHS FAS HM RD FADER AH HY ERK HH EH GEE HATE AEA RK ESE GS ERE ERS AE THERES HED 


€BLCCK 8 CONTRCL PARAMETERS 3 CAROS FORMAT a110 
: NEB ; NEF NJ ; NCJ NMT TOVCLC NOJ 
2 | 2 5 3 2 2 S 
: NEUBC LBLCK l NFREQ NFMASS , NEIG l NEIGL NPRI 
il 1 1 
: NLC ' NCSFLC ' NSTRES ; NSTW NF MW 
2 
eres iS ie a ss Paty 50d Sieh ERS S$ 0Gk 24 052,406 EnSeS 04 6 XERKG SES Adasen SETS 784 EO SRO 
$2LOCK C DYNAMIC ANALYSIS INFORMATION 1 CARD FORMAT 110,2F10.0 
§ LMASS GREV SPSeIG 
- Q : : 
LR SREKB ID] KM BY eK OEE EEK SZC HAAEREHRGE SH SEEK ARES EHR SEKESE HAK GHAR KAT RKATHHAEE 
f@tCCK C JOINT CCORCINATES 
s NJ CAPOS FORMAT I[10-¢3F10 
: , X-CCCRO Y-COCRO Z-C CORO 
1 Oo O~ Oe 
2 100. orn Oe 
3 2006 Oe O~« 
4 Oo 150. Oe 
5 Oo 0. 50. 
$8LCCK £& DESIGN VARIABLE LINKING DATA 
¢ mJ) CAROS FORMAT 4110,3F10.0 
$ JN [x ly [TZ PCX PCY PCZ 
$ { | { | { { 
~ C Q 120 1.9 1.0 
5 0 8) 1.0 1.0 1.0 
$ SBR BRAY] ERK AVS HK AD EK ASHKHE SH RAK HK KEG ASE ERAS RARE ASEK ERA RARER] RRKEKKSSSE 
tB8LOCK F MATERTAL PRCPERTIES NMT CAROS FORMAT 6F10290 
2 = REC SIGMIN SIGMAX K EULER POISSN 
eae OO. -25000. 25CO00.% 24 oh 
Ze SE+7 0.3 -36000. 36000. ° Beat 
EASRAKRA BAL SRA AVHLRM RIRKK AK SAKE A SHEAKA TEAGAT HOSA AAR OME RASRA SAA SAAALAASS 7 KABA ORE 
$BLCCK G GAR ELEMENT INFORMATION NEB CAROS FORMAT 5110¢F10,110 
: LNO NCCE2 NCDE2 MATCOC NSO0G1 | AREA | ESeCl 
1 b - 1 1 2.0 1 
2 2 5 : 2 2.0 1 
SSSSREAGA]L SAARAGSSS SD SAPVSAKRSASZ SS SHKARARER4GKS FR SK SKS KIER KARL EKETRASE ARASH TAKERS FEB 
SSPeCK snore AM: ELEMERT IN-e-CGrRMA: SON Z2*NEF CAROS FORMAT 71I10/2F 10 
: LNO M@IGE 2 NCDE3 MATCOD NSDG1 NSOG2 LSecT 
S 2 2 %) + 
i ees 
£.0 1.9 
4 2 3 Pa a 5 di 
ee 8 Neo Yk. ae an Send ae de Gearie dd S¥2cb S82 242 x6 2 ee ee seat bees 
seueeck I JCINY CGNSTRAINT OATA NCJ CAROS FORMAT 7110 7 
: JN [x ry IZ IXX Leyy. bee 
L { 
2 0) 8) 0 0 0 0 
3 C 0 0 8) e) 0 
4 1 l i 1 - ; 
Deere F0ys (54 5 aon 24.99 4g 20s 944-6 Se EEE Geo 2 oign We bese a 940 54 5e ONES RT +2 Fees A e8 
SelGCK J JCINT LGACING CATA 1 CARO FORMAT 110 
$ GGyits:f NUC=C IN BLOCK 8B) 
$ NLJ 
$ { 
1L 





TABLE “XT 


DATA FILE TRUSS-BRACED CANTILEVER BEAM continued 


$ NLJ CAROS  FORMATI10,6F10 
a. FX i / FL TK a Tm 
3 0. 1000.0 de 0. O. or 
il 
3 0. 0.0 -1000.0 0. 0. 0. 
ge oeses ea] te esaeta st 244954558 x2 eet aes NEA eRe oA See ed Seen nt caer aS 
SBLGCK K FIXED MASS CATA NEMASS CARDS EGRMAT 110,F10.0 
g (CMIT TF NEMASS20 IN BLOCK B) 
$ JN MASS 
/ h 250.4 
SSRIS RH AM] SRR ADH II SRE RSH AAKSEHSE aK EKXRGHK RRSKEKKSHS MR AREKKKRE EE SEAR AAKEKKI ORES kK es 
SBLCCK L DESIGA VARTARLE INFORMATION CaRDS AS REQ!) FORMAT 8F10.0 
$ (AREA ANC CIMENSISN VAPIABLES CMIT IF NOVAP1=0) 
$ XACT) «a DA CNOVAR | | , , 
2C6C 20-0 200 2.0 2.0 
$  XAL(T) 
$ | { ! | | j { 
250 - 50 400 210 210 
$ xaull) 
$ | | { | | j 
2£.0 aGec 25.0 2.5 4.0 
tH AVaaRAS)] SaHS BRR EzSeexeus ABE SX Sex 4xe AK AES RAS EEEREOSAEEEE RES TX EES 
SBLCCK M DESIGA VARIABLE INFORMATION CARDS AS REQ'D FORMAT 8F1020 
f (CCORDINATE VARIGELES SME? 12 NBVAR2=0) 
$ CUT} ee XC{NGVAR 2 } 
¢ i | | | ' 
150.0 £0.20 
$  xCL{t) 
$ | | } j \ j | 
A xcuis Ero 
- 
200.0 100 .0 
x SREKKRL] XEXS RAMEE SRE SMR ee GOK EK FAAS RAKE EMAAR ASS eee TES EES x3e 
€BLCCK N JCINT CISPLACEMENT CONSTRAINTS NDSPLC CARDS FCRMAT 3110,2F10 
g (CMIT IF KOSPUC=c IN BLGCK 3B) 
¢ JN CIR LC BU 
$ | ) | | | 
2 5 = 73-9 zee 
CeRSSe DHT Se RSVR %*x2daua%4~% 33 aes ERX ase, ke xx xe we KS “xx ARR KS See TeHEKERSY za 
SBLCCK C FREQUENCY CCNST2AINTS KO. OF CARDS AS REQ'D ECRMAT 8F10.0 
‘ (CMIT TF NFREG=0 IN BLOCK 3} 
$ FC s @ C} 
¢ 7. | i j i | | 
TTTT TESTE TICS TPL ITIL COEcCe CT Te et Teer ere Te Tie TE TCCCSTCe TT reres + Pe Te STL eee 
SELCCK F ENO CARD 


~] 
OO 





TRB We 
ween Ley TWO}TLER 3=D “PORPAL FRAME 


SSK SHEEHY) SE TIDES HS EEE SZ BREE EK EKG OHS ARS ERASE EL EHREHERAEEET EEREEE NS 
marek a me 1 CARO FCRMAT 80H 
-~TIEF FRAME STRUCTURE WITH TRUSS X=8RACES: EXAMPLE 2 


ee SOG ees Co So ae chee aes Bene se BAER SEARS REE AKEEAGRAAEKKAARK] AKKGSKSE EKS 
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PROGRAM ORGANIZATION 


A. DESCRIPTION 

Brew progrem organizat#on is layed out in the following 
FYow Charts. The main driver program (SADXM) calls a 
subdriver (SADXSD), and the optimizer of the user's choice. 
AVI Changes required for replacement of the optimizer are 


Ma@e in SADXM, thus allowing for easy testing of several 


PEI zers on the same problem. 


Sense mey be called from the main @or input, analysis, 
Pea OUspe.. Brinked output may vary as the user requires. 
A complete listing of all subroutines and their functions is 
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